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The composite rotation approach has been used to develop a range of robust quantum logic gates, 
including single qubit gates and two qubit gates, which are resistant to systematic errors in their 
implementation. Single qubit gates based on the BBl family of composite rotations have been exper- 
imentally demonstrated in a variety of systems, but little study has been made of their application in 
extended computations, and there has been no experimental study of the corresponding robust two 
qubit gates to date. Here we describe an application of robust gates to Nuclear Magnetic Resonance 
(NMR) studies of approximate quantum counting. We find that the BBl family of robust gates is 
indeed useful, but that the related NBl, PBl, B4 and P4 families of tailored logic gates are less 
useful than initially expected. 
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Quantum information processing [l| has made substan- 
tial progress in recent years, but formidable problems 
remain in the implementation of a general purpose quan- 
tum cornputer. The development of quantum error cor- 
rection Q, d, Q was a key step, as it allows quantum 
computers to function in the presence of random errors. 
More recently the method of composite rotations (also 
called composite pulses) , developed for Nuclear Magnetic 
Resonance (NMR) experiments ^i^iTj,^), has been used 
to design quantum logic gates which are robust against 
systematic errors in the control fields used to implement 
them ^iia]. 

Composite pulses are closely related to the more gen- 
eral field of optimal quantum control, which has been 
widely studied and finds applications in many fields. In 
the context of conventional NMR experiments there has 
been much study of the design of time-optimal pulse se- 
quences , which enable unitary transformations to be 
performed while minimizing losses from dccoherence. For 
applications in NMR quantum computing there has been 
particular interest in the use of strongly modulated com- 
posite pulses 0, 0, 01 , which are designed to perform 
particular selective operations in the presence of a com- 
plex multi-spin Hamiltonian, and which can also be de- 
signed to be robust against systematic errors in control 
fields. Although often discussed separately there are clear 
similarities between these approaches |l5j . Similar meth- 
ods have also been developed in ion trap implementations 
of quantum computing 16]. 

Henceforth we only consider a particularly simple 
group of composite gates, based on the composite pulses 
originally developed for conventional NMR experiments. 
These pulses differ from the more general approaches 
mentioned above in two significant ways. Firstly they 
provide a simple analytic set of solutions, which allow 
simple rotations to be converted directly into composite 
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rotations, rather than a recipe for finding a numeric solu- 
tion for a particular problem. Secondly, the implementa- 
tion of these composite gates is particularly simple. For 
example, in NMR systems they only require control of the 
phase and duration of a small number of RF pulses at a 
single fixed frequency and power, while strongly mod- 
ulated composite pulses also require the frequency and 
power to be varied. For this reason, these pulses are easy 
to implement experimentally. 

Single qubit gates developed using this approach have 
been demonstrated experimentally in a range of systems, 
including NMR Electron Spin Resonance (ESR) 

[Ulill, and the Quantronium SQUID but there 

has not yet been much study of their use in an extended 
quantum algorithm. A related family of robust two qubit 
gates has also been described , but these partic- 

ular gates have not so far been studied experimentally. 
Here we described the application of single and two qubit 
logic gates based on the BBl, NBl and PBl families of 
composite pulses to an implementation of approximate 
quantum counting using NMR studies of a heteronuclear 
spin system. 



I. THE BBl, NBl AND PBl FAMILIES 

Composite pulses are designed to perform quantum op- 
erations in the presence of systematic errors. Note that it 
is not necessary to know the size of the error involved, but 
it is necessary to know its general form. They have been 
extensively developed in conventional NMR studies 
principally to correct pulse length errors, that is errors in 
the strength of the field used to induce a rotation, and 
off-resonance errors, that is imperfections arising when 
an oscillating field is not quite in resonance with the tran- 
sition it is supposed to drive. Most of these composite 
pulses are not suitable for use in quantum information 
processing, but a small group of them, sometimes called 
fully compensating pulses, are well suited. These pulses 
work on any initial state, and can in principle be used to 
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replace naive gates without any further modifications. 

Perhaps the most useful family of fully compensating 
pulses developed to date is the BBl family of pulses which 
are robust against pulse length errors. The ideal unitary 
transformation for a 0^ pulse is 

U{9, 0) = exp[-i0{lx cos + Jy sin 0)] (1) 

but in a real physical system the actual operation will not 
have this perfect form. In the presence of pulse length 
errors the real operation has the form 

¥{9, 0) = exp[-i{l + f)e{I, cos c^ + Iy sin cj))] (2) 

where / is the fractional error in the pulse length or field 
power. Defining the propagator fidelity as 



\Tr{VU^)\/Tr{UU^) 



gives the result 



i^ = cos(^) 
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for a simple 6*^ pulse. The BBl family of robust pulses 
is obtained by replacing the simple pulse 9^ by the com- 
posite pulse sequence 

{9/2)^{tt)^+^{2tt)^+3^{tt)^+^{9/2)^ (5) 

where ip — arccos(— 6'/47r). This composite pulse has a 
fidelity of 



6 {3271^9^ + U7T^9^ - 9^) 
9216 



(6) 



showing that for small errors the BBl pulse is far less 
error prone than a simple pulse. 

The NBl family of pulses is obtained in a similar way; 
in this case the simple pulse is replaced by 



(6l/2)0(7r)0+^(27r)0_^(7r)0+^(6'/2)0 



(7) 



which has the same structure as BBl, differing only in the 
phases. This sequence is more error-prone than a simple 
6*0 pulse, and so cannot be used to suppress errors. It 
docs, however, have the interesting property that when 
I/I « 1 the resulting composite pulse is a good approxi- 
mation to the Identity operation. Thus is the case of very 
low pulse powers (/ w —1) an NBl pulse "does nothing" : 
in effect, evolution under weak fields is suppressed. 

The PBl family of pulses offers a compromise between 
these two extremes. The simple pulse is replaced by 



{9/2)^{2n)^+^, (47r)^_^, {2tt)^+^, {9/2)^ 



(8) 



where ip' = arccos(— ^/Stt). These pulses are both ro- 
bust to small errors and able to suppress weak fields, 
suggesting that PBl is the best general purpose family 
of pulses, although the BBl and NBl families perform 
their respective tasks more effectively. 



Since these families were described, a range of similar 
pulses have been discovered, most notably the arbitrary 
precision pulses of Brown et al. jU|2^. In their notation 
BBl is called B2 and PBl is called P2; these are the 
first members of a series of families of pulses with ever 
greater error tolerance. These sequences swiftly become 
extremely long, and here we consider only the next two 
members, B4 and P4. The B4 composite pulse takes the 
form 



{9/2)^ [(7r)0+^(27r)0+3^(7r)0+^]'' 

(— 27r)0+^(— 47r)0_^,(— 27r)0+^ 

[(7r)0+i>(27r)0+3^(7r)0+^]** (9/2)^ (9) 

where ip = arccos(— 0/247r) and the superscript 4 indi- 
cates that the section enclosed by square brackets is re- 
peated four times. The central three pulses have negative 
rotation angles; these can be partially canceled with sur- 
rounding pulses, and the remaining rotations can be im- 
plemented as positive rotations around axes with phase 
angles offset by tt. The P4 composite pulse is similar, 
taking the form 

(9/2)^ [i2n)^+^, (47r)^_^. i2n)^+^,t 

(— 47r)0+^'(— 87r)0_^'(— 47r)0+,0/ 

[(2^)^+0, (4^)0_^,(27r)^+0,]* (9/2)^ (10) 



where ■tp' = arccos(— 6'/487r). Note that the forms for 
these pulses origii iallv published by Brown et al. |22| are 
slightly incorrect 

Under ideal conditions B4 and P4 perform slightly bet- 
ter than their simpler counterparts, but the gain is slight 
and more detailed simulations [2J| indicate that these 
pulses are highly sensitive to the presence of off-resonance 
and phase errors. This is confirmed by our experience, 
described below, that the performance of B4 and P4 is in 
practice quite poor. Thus the original three families offer 
perhaps the best combination of simplicity and effective- 
ness, and we concentrate on them in most of what follows. 
It might seem that, as previously suggested, PBl is the 
best general purpose family of pulses, but this turns out 
not to be the case. 



II. THE QUANTUM COUNTING 
EXPERIMENT 

Most studies of these particular robust quantum logic 
gates to date have been either theoretical or have in- 
volved only simple demonstrations of single logic gates. 
Although these are of some interest, it is also important 
to investigate more complex situations, such as quan- 
tum algorithms containing large numbers of logic gates. 
Quantum counting provides an ideal testing ground as it 
permits large numbers of logic gates to be implemented 
in a simple physical system. 

Quantum counting is closely related to Grover's quan- 
tum search [25l l2a . l27 | . Consider a function f{x) which 
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FIG. 1: A quantum circuit for implementing quantum count- 
ing on a two qubit NMR quantum computer; the central se- 
quence of gates, surrounded by brackets, is applied r times. 
A similar circuit can be constructed for a larger search space 
by replacing the (lower) target bit by a register and replac- 
ing gates applied to the target by multi-bit versions. Gates 
marked h implement the NMR pseudo-Hadamard operation, 
while those marked implement the inverse operation. 




maps n-bit strings to a single output bit, such that 
f{x) = or 1. In general there are N — 2" possible 
input values, with k values for which f{x) — 1. Grover's 
quantum search allows one of these k items to be located, 
while quantum counting |2p| allows the value of k 
to be estimated. The counting algorithm estimates an 
eigenvalue of the Grover iterate G = HUqH^^Uj, which 
forms the basis of the searching algorithm, where H is 
the rt-bit Hadamard transform, Uq maps |000 ... 0) to 
— 1000 . . .0) and leaves the remaining basis states alone, 
and C/jmaps to (— For further details see 

A quantum circuit which implements this algorithm on 
a two qubit NMR quantum computer is shown in Fig. ^ 
As usual, pairs of Hadamard gates have been replaced 
by NMR pseudo-Hadamard gates and their inverse pol 
|30j. This circuit can be used to count the number of 
solutions to f{x) = 1 over a one bit search space, but 
similar circuits exist for larger search spaces. 



III. SINGLE QUBIT GATES 

To investigate the effects of systematic errors on sin- 
gle qubit gates we implemented this circuit in an NMR 
experiment using a solution of ^'^C labeled sodium for- 
mate in D2O [31|. The formate anion HCO^ provides 
an isolated two spin system, comprising a ^H and a ^^C 
nucleus, with a large coupling between the nuclei. This 
coupling arises from the scalar coupling interaction, and 
in a heteronuclear spin system takes the weak-coupling 
form, usually written as 7rJ2IzSz in NMR product op- 
erator notation Q. (This coupUng is sometimes said to 
have the Ising form, although this name more properly 
refers to extended networks of such couplings.) In a het- 
eronuclear spin system the RF transmitters can be place 
in perfect resonance with each qubit, and so off-resonance 
effects can be essentially ignored. The major remaining 
source of systematic error is pulse length errors arising 
from inhomogeneity of the RF fields. 



FIG. 2: The effects of pulse length errors on a single qubit 
pseudo-Hadamard gate (a 90^ rotation). The figures show the 
effect of applying this gate to a spin in the thermal equilibrium 
state, which is equivalent to a pseudo-pure state |0), followed 
by observation of the NMR spectrum, effectively measuring 
the off-diagonal components of the single-spin density matrix. 
For further details see the main text. 



We begin by demonstrating the effect of using robust 
quantum logic gates to implement a pseudo-Hadamard 
gate (a 90° rotation) on the ^H qubit which begins in 
an equilibrium state. In addition to the pulse length er- 
rors which arise naturally from inhomogeneity, additional 
errors can be introduced by deliberately mis-setting the 
pulse length. The effect of doing this is shown in Fig. |2 
Each sub-figure contains a set of spectra corresponding 
to pulse length errors in the range ±100%, while the six 
sub-figures correspond to naive pulses and the five fami- 
lies of composite pulses described previously. 

As expected the naive pulse gives a signal whose inten- 
sity is a cosine function of the fractional error. The BBl 
family is much more robust to pulse length errors, while 
the NBl family is more sensitive to errors than the naive 
pulse, and gives much smaller excitations than the naive 
pulse for errors close to ±100%. The PBl family shows 
the expected compromise behaviour, with a broad cen- 
tral maximum for small errors, and broad minima around 
errors of ±100%. Note that the performance is better for 
errors around —100%, which is the experimentally impor- 
tant case of very weak fields, than for the theoretically 
equivalent but experimentally unimportant case of errors 
close to +100%; this reflects the intrinsic inhomogeneity 
of the RF field. The performance of the B4 composite 
pulse is broadly similar to that of the much simpler BBl 
sequence, while the performance of P4 is clearly rather 
poor. 

This shows that the robust logic gates can have the 
desired effects when used to implement a single logic 
gate, with BBl and PBl offering the best performance, 
but it is also important to investigate how they work in 
more complex situations. This was done by implement- 
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FIG. 3: An implementation of approximate quantum counting with naive and BBl single qubit gates. The use of BBl single 
qubit gates greatly reduces the apparent decoherence rate, indicating that much of the signal loss actually arises from pulse 
length errors, and that BBl is effective in correcting for this in a complex sequence. 



ing the quantum counting circuit with the results shown 
in Fig. 121 Only naive and BBl single qubit gates were 
used. The improvement obtained by using BBl single 
qubit gates is clear, with much of the signal loss that 
would naively be ascribed to decoherence clearly arising 
from pulse length errors. BBl single qubit gates are, 
therefore, used throughout the remainder of this paper. 



IV. TWO QUBIT GATES 

A very similar approach can also be used to tackle sys- 
tematic errors in coupling gates, which provide the ba- 
sic two qubit gate for NMR quantum computation |30j |. 
Evolution under a scalar coupling can be thought of as a 
rotation about the 2IzSz axis, and errors in the coupling 
constant J correspond to errors in a rotation angle about 
this axis. Such errors can be overcome [l^,!^ by rotat- 
ing about a sequence of axes tilted from 2IzSz towards 
another axis, such as 2IzSx- Defining 

6*0 = exp[~i X0X {2IzSz cos0 + 2IzSx sin0)] (11) 

allows the naive sequence to be replaced by Eq. as 
before. The tilted evolutions can be achieved by sand- 
wiching free evolution under the coupling Hamiltonian 
between (f>z^y pulses applied to spin S. For the case that 
6 ^ tt/2 (which forms the basis of the controlled-NOT 
gate) the final sequence takes the form shown in Fig.^ 

The results of using these robust two qubit gates could 
be investigated in much the same way as the single qubit 
gates shown above, but it is more interesting to con- 
sider another phenomenon, that is the ability of the NBl 
and PBl gates to suppress evolution under small scalar 
couplings. This was studied by implementing quantum 
counting on the more complex spin system provided by 
alanine which is ^'^C labeled at position 2. Like the for- 
mate anion, this contains a CH spin system, comprising 
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FIG. 4: Pulse sequence for a BBl robust coupling gate to 
implement a controlled-NOT gate in a two spin (IS) system. 
Boxes correspond to single qubit rotations with rotation an- 
gles oftp — arccos(— 1/8) ~ 97.2° applied along the ±j/ axes as 
indicated; time periods correspond to free evolution under the 
scalar coupling, nJ2IzSz, for multiples of the time t = 1/4 J. 
The naive coupling gate corresponds to free evolution for a 
time 2t. 



the labeled ^^C and the directly bonded nucleus, but 
unlike formate this spin system is not completely isolated, 
as both nuclei have small couplings to the nuclei in 
the methyl group. This additional coupling interaction 
could be removed by deuterating the methyl group, or by 
applying decoupling fields, but here we seek to suppress 
the coupling using robust logic gates. 

We begin by studying evolution under the CH cou- 
pling, and concentrate on evolution times of n/2J, where 
n is a positive integer. This was done by applying a 
pseudo-Hadamard gate to the ^'^C nucleus and then al- 
lowing the spin system to evolve under a coupling gate 
(either naive or composite) for an effective time n/2J. 
The results are shown in Fig. [S] The naive gate appears 
to perform well for small evolution times, as the evolu- 
tion under the small J coupling can be largely ignored: 
this coupling is clearly visible in the spectra where each 
component of the doublet is split into a 1 : 3 : 3 : 1 multi- 
plet, but all four components appear to have very similar 
phases. For longer evolution times, however, the detri- 
mental effect of this additional coupling becomes clear, 
as the components appear with quite different phases, 
indicating significant evolution under the coupling. 
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FIG. 5: Implementation of an coupling gate in the presence of extraneous couplings. For details see the main text. The ten 
^^C spectra in each subfigure correspond to increasing periods of evolution under the coupling, and the desired pattern is a 
simple alternation in signal intensity between +1 and —1; deviations in the naive implementation arise from evolution under the 
extraneous couplings to the methyl protons. It might be expected that evolution under this small coupling would be effectively 
suppressed by the NBl sequence, but this performs very poorly, and the best results are seen from the BBl sequence. 



In an attempt to overcome this we implemented NBl 
coupling gates, which should suppress evolution under 
this small coupling. It is immediately obvious from Fig.jS] 
that this approach does not work, as the multiplets are 
far more distorted than those seen using the naive se- 
quence. In retrospect this is unsurprising: NBl is de- 
signed to suppress a small coupling on the assumption 
that this is the only coupling present. If the small cou- 
pling occurs in addition to another much larger coupling 
then there is no reason to believe that NBl will suppress 
it. 



Examination of the results from BBl and PBl cou- 
pling gates makes this point even more clearly. BBl gives 
much cleaner spectra than the naive coupling gate, with 
the results of PBl being similar but significantly poorer. 
Once again the explanation is clear in retrospect: al- 
though we know that the splittings visible in the spec- 
trum arise from the combination of a large coupling and 
a small coupling, the pulse sequence is blind to this ori- 
gin. The same pattern could in principle arise from a 
mixture of different molecules with a range of coupling 
constants. The BBl robust gate is designed to give very 
similar behaviour over this range of coupling sizes, and 
so all the different components of the multiplet appear 
inphase. The action of PBl is similar: the quality of the 
spectra seen arises not from the ability of this sequence 
to suppress small couplings, but from its ability to tol- 
erate a range of couplings. Over the range of couplings 
seen in this system the behaviour of PBl and BBl should 
be very similar, and the relatively poor results observed 
from PBl are probably a consequence of the fact that the 
PBl sequence take almost twice as long as BBl, leading 
to increased problems from decoherence. 



V. SIMPLIFYING SPIN SYSTEMS 

The results above show that our original idea that NBl 
based composite two qubit gates could be used to sup- 
press small couplings, in effect simplifying complex spin 
systems, is incorrect. However, BBl based two qubit 
gates do result in the desired suppression, with PBl se- 
quences giving similar but slightly poorer results, sug- 
gesting that it might be possible to use these to simplify 
spin systems instead. 

Following this idea we attempted to implement quan- 
tum counting in our labeled alanine spin system, using 
BBl to suppress the small couplings between the ^^C 
qubit and the methyl protons while using the large cou- 
pling between the ^■^C and the directly bonded qubit 
to implement logic gates. These attempts were not suc- 
cessful (data not shown), and once again it is clear in 
retrospect why this idea will not work. Although the 
small couplings between the ^^C and the methyl pro- 
tons in the labeled alanine spin system can indeed be 
suppressed by BBl, there are also small couplings be- 
tween the qubit and the methyl protons which will 
not be suppressed. These homonuclear couplings could, 
of course, be removed by frequency selective pulses, or 
by decoupling, but either approach would act to directly 
simplify the spin system, rendering the BBl approach 
unnecessary. 



VI. CONCLUSIONS 

Our results confirm that simple composite pulses can 
indeed be used to suppress systematic errors on single 
qubit gates used in implementations of complex quantum 
algorithms. The BBl approach is likely to prove the most 
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useful, while PBl may find applications in some special 
areas. The BBl sequence can also be used to implement 
robust two qubit gates, although in this case the extra 
time required for the longer pulse sequence may cause 
difficulties. The more complex B4 and P4 sequences, 
although theoretically superior, do not perform well in 
practice. 

Two qubit gates based on NBl and PBl could in prin- 
ciple be used to suppress the effects of small couplings, 
but this is not effective when the small couplings occur in 
addition to larger couplings which are used to implement 
gates. In this case the small couplings can be treated as 
small errors in the large coupling, and BBl provides the 
best suppression of their effects. It might appear that 
this approach could be used to neglect small couplings 
in complex spin systems, effectively simplifying the spin 



system through composite pulses, but this approach will 
rarely if ever be effective, as corresponding homonuclear 
couplings cannot be suppressed by these simple compos- 
ite pulse methods. More complex approaches, such as 
strongly modulated composite pulses, could be more ef- 
fective, but such pulses need to be designed on a case by 
case basis. 
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